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Abstract 

The algebra of fourvectors is described. The fourvectors are more appropriate than the 
Hamilton quaternions for its use in Physics and the sciences in general. The fourvectors 
embrace the 3D vectors in a natural form. It is shown the excellent ability to perform 
rotations with the use of fourvectors, as well as their use in relativity for producing 
Lorentz boosts, which are understood as simple rotations. 
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1 Introduction 

1.1 General 

In this paper it is suggested the use of 
fourvectors with the purpose of replacing the 
3D vectors and the quaternions. Because 
the fourvectors contain the three dimensional 
vectors and can be considered a formalization 
of them. 

The discovery of the quaternions is attributed 
to the Irish mathematician William Rowan 
Hamilton in 1843 and they have been used for 
the study of several areas of Physics, such as 
mechanics, electromagnetism, rotations and 
relativity [26], [20], [6], [2], [23], [13]. James 
Clerk Maxwell used the quaternion calculus 
in his Treatise on Electricity and Magnetism, 
published in 1873 [22]- An extensive bibliog- 
raphy of more than one thousand references 
about Quaternions in mathematical physics 
has been compiled by Gsponer and Hurni 

The modern vectors were discovered by 
the Americans Gibbs and Heaviside between 
1888 and 1894. Their work may be consid- 
ered a sort of combination of quaternions and 
ideas developed around 1840 by the German 
Hermann Grassman. The notation was pri- 
marily borrowed from quaternions but the 
geometric interpretation was borrowed from 
Grassman's system. 

By the end of the nineteenth century 
the mathematicians and physicists were hav- 
ing difficulty in applying the quaternions to 
Physics. 

Ryan J. Wisnesky [32] explains that "The 
difficulty was a purely pragmatic one, which 



Heaviside was expressing when he wrote that 
there is much more thinking to be done [to set 
up quaternion equations]. In principle, most 
everything done with the new system of vec- 
tors could be done with quaternions, but the 
operations required to make the quaternions 
behave like vectors added difficulty to using 
them and provided little benefit to the physi- 
cist." 

"Gibbs was acutely aware that quater- 
nionic methods contained the most important 
pieces of his vector methods." [32] 

After a heated debate, "by 1894 the debate 
had largely been settled in favor of modern 
vectors" [22]. 

Alexander MacFarlane was one of the 
debaters and seems to have been one of the 
few in realizing what the real problem with 
the quaternions was. "MacFarlane's attitude 
was intermediate - between the position of 
the defenders of the Gibbs-Heaviside system 
and that of the quaternionists. He supported 
the use of the complete quaternionic product 
of two vectors, but he accepted that the 
scalar part of this product should have a 
positive sign. According to MacFarlane the 
equation j k = i was a convention that 
should be interpreted in a geometrical way, 
but he did not accept that it implied the 
negative sign of the scalar product". [25] 
(The emphases are mine). 
He incorrectly attributed the problem to a 
secondary and superficial matter of repre- 
sentation of symbols, instead of blaming to 
the more profound definition of the quater- 
nion product. "MacFarlane credited the 
controversy concerning the sign of the scalar 
product to the conceptual mixture done by 
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Hamilton and Tait. He made clear that the 
negative sign came from the use of the same 
symbol to represent both a quadrantal versor 
and a unitary vector. His view was that 
different symbols should be used to represent 
those different entities." [25] (The emphasis 
is mine). 

At the beginning of the twenty century, 
Physics in general, and relativity theory 
in particular, was lacking an appropriate 
mathematical formalism to represent the 
new physical quantities that were being 
discovered. But, despite the fact that all 
physical variables such as space-time points, 
velocities, potentials, currents, etc., were rec- 
ognized that must be represented with four 
values, the quaternions were not being used 
to represent and manipulate them. It was 
necessary to develop some new mathematical 
devices to manipulate such variables. Be- 
sides vectors, other systems such as tensors, 
spinors, matrices and geometric algebra were 
developed or used to handle the physical 
variables. 

During the twenty century we have wit- 
nessed further efforts to overcome the diffi- 
culties remaining, with the development of 
other algebras, which recast several of the 
ideas of Grassman, Hamilton and Clifford in 
a slightly different framework. An example in 
this direction is Hestenes' Geometric Algebra 
in three dimensions and Space Time Algebra 
in four dimensions. [16], [T7], [21], [8] [18] 

The commutativity of the product was 
abandoned in all the previous quaternions 
and in some algebras, such as the one of Clif- 



ford. According to Gaston Casanova [I] "It 
was the English Clifford who carried out the 
decisive path of abandoning all the commu- 
tativity for the vectors but conserving their 
associativity." "This algebra absorbs the 
Hamilton quaternions, the Girard's complex 
quaternions, the cross product and the com- 
plex numbers, the hyperbolic numbers and 
the dual numbers." [4]. Also the Hestenes' 
"geometric product" conserves associativity 
[IS] . In this sense, the associativity of the 
product is finally abandoned in the fourvector 
algebra proposed in the present paper. This 
means that the fourvectors do not constitute 
a Clifford Algebra [2] or a Geometric Alge- 
bra pp. This is a collateral effect of the pro- 
posed algebra, and constitutes a hint about 
the form the fourvectors should handle, for 
example, a sequence of rotations; besides, the 
complex numbers are not handled as in the 
Hamilton quaternions, where the real number 
is put in the scalar part and the imaginary in 
the vector part, but a whole complex number 
can be put in each component, so it is pos- 
sible to have up to four complex numbers in 
each fourvector. But, what is more impor- 
tant, it is known that in quantum mechanics, 
observables do not form an associative alge- 
bra, so this could be the natural algebra for 
Physics. 

The proposed algebra could have been 
already developed, around 1900, under 
the name of hyperbolic quaternion, which 
is a mathematical concept introduced by 
Alexander MacFarlane of Chatham, Ontario. 
The idea was dismissed for its failure to 
conform to associativity of multiplication, 
but it has a legacy in Minkowski space and 
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as an extension of "split-complex numbers". 
Like the quaternions, it is a vector space 
over the real numbers of dimension 4. There 
is only the mention to such quaternions but 
no accessible references to confirm if those 
quaternions satisfy the same algebraic rules 
given in the following for the fourvectors. 
However our intent is to convince the reader 
that the presented here is one of the most 
important mathematical tools for Physics. 

The fourvector representation is without a 
doubt a more unified theory in comparison to 
the classical vector or matrix representations. 

The use of fourvectors allows discerning 
constants, variables and relations, previously 
unknown to Physics, which are needed to 
complete and make coherent the theory. 

The vectors have lost some ground in 
favor of the Hamilton quaternions due to 
the lack of an appropriate 4D-algebra. For 
example, Douglas Sweetser, who has worked 
extensively in the application of Hamilton 
quaternions to many possible physical areas, 
in general with very little success, sustains 
these opinions: "Today, quaternions are of 
interest to historians of mathematics. Vector 
analysis performs the daily mathematical 
routine that could also be done with quater- 
nions. I personally think that there may be 
4D roads in physics that can be efficiently 
traveled only by quaternions." [29] 

In fact those 4D roads should be traveled 
only by properly handled fourvectors. It has 
been an old dream to express the laws of 



Physics using quaternions. But this attempt 
has been plagued with recurring pitfalls for 
reasons until now unknown to both physi- 
cists and mathematicians. The quaternions 
have not been making problem solving easier 
or simplifying the equations. Very often the 
Hamilton quaternions require an extreme lia- 
bility to guess when and where a quaternion 
needs to be conjugated, in order to obtain 
some particular result. 

I believe that this has been due to an 
internal problem in the mathematical struc- 
ture of the Hamilton quaternions, which I 
will try to reveal in the present paper. The 
correction of such problem constitutes a new 
non-commutative, non-associative normed 
algebraic structure with which it is possible 
to work with fourvectors in an improved 
way relative to the Hamilton, Pauli or Dirac 
quaternions, geometric algebra, space-time 
algebra and other formalisms. 

In the present paper, in particular, the 
present author exhibits the application of the 
fourvectors to 3D and 4D rotations, which 
requires a reformulation of the Hamiltonian 
mathematics. 

The powerful Mathematical package 
includes a standard algebra package for the 
manipulation of the Hamilton quaternions. 
I have borrowed from that package the 
symbol, as double asterisk, to represent the 
fourvector product. It is easy to modify the 
cited package to handle the fourvectors, as 
well as to permit not only their numeric but 
also symbolic and complex manipulation. 
Though I have still not been able to figure 
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out a simple way to reproduce the trigono- 
metric and exponential functions for the 
fourvectors (if at all possible), which that 
package allows to compute for the Hamilton 
quaternions. 

In the following three subsections, a cur- 
sory revision is made of the Hamilton, Pauli 
and Dirac quaternions, for an easy compari- 
son with fourvectors. In section 2 the fourvec- 
tors are presented. Section 3 can be skimmed 
by the mathematician, since it is mostly clas- 
sic algebra. Finally, in section 4 the formulae 
needed to perform rotations and reflections 
with fourvectors is presented. 

1.2 The Hamilton quaternions 

Quaternions are "four-dimensional numbers" 
of the form [31 J: 

A = a + \a x + j a y + ka 2 , (1) 
B = b + i b x + j by + k b z 

where the basis elements 1, i, j, k satisfy the 
relations: 

i 2 =j 2 = k 2 = ijk=-l (2) 
and also: 

ij = -ji = k, (3) 

jk=-kj = i, (4) 

ki = -ik = j. (5) 

Here 1 is the usual real unit; its product 
with i, j or k leaves them unchanged. 
Thus, since the products of the basis elements 



are non-commutative, we have in general 
A**B ^ B**A, where the double aster- 
isk represents quaternion multiplication. Un- 
der these conditions, quaternion multiplica- 
tion is associative, so that (A**B)**C = 
A**(B**C) for any three quaternions A, B, 
C. 

The sum of two quaternions is 
A + B =(a + b) +i(a x + b x ) + 

(b) 

j(a„ + by) + k(a 2 + b z ), 

and, using relations (j2J) and (jSI)-©, the prod- 
uct is given by: 

A * * B = {ab — a x b x — a y b y — a z b z ) + 
i (ab x + a x b + a y b z - a z b y ) + 
j (ab y - a x b z + a y b + a z b x ) + 
k (ab z + a x b y - a y b x + a z b). 

The notation of three-dimensional vec- 
tor analysis furnish a useful shorthand for 
quaternion operations. Regarding i, j, k as 
unit vectors in a Cartesian coordinate system, 
we interpret the quaternion A as comprising 
the scalar part a and the vector part 
a = ia I +ja l/ + ka 2 . Then we write in the 
simplified form A = (a, a). With this nota- 
tion, the sum (j5J) and the product (j7j) may 
more compactly be expressed as: 

A + B =(o+ fe,a + b) (8) 

A**B=(a b — a ■ b, ab + 6a + ax b) 

(9) 
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where the usual rules for vector sum and 
dot and cross products are being invoked. 

According to the mathematicians, the 
Hamilton quaternions are mathematical 
structures which combine properties of com- 
plex numbers and vectors. 

1.3 The Pauli quaternions 

The Hamilton multiplication rules differ from 
the Pauli matrices rules only by the explicit 
appearance of the fourth basis element. 

The basis elements of the Pauli quaternion 
space are denoted by Si, s 2 , s 3 , s 4 . 

They obey the following multiplication 
rules, comparable to (0)-©: 



2 2 
S l — S 2 


= 4 = 




s l s 2 — 


— S 2 Si 


= S3 


S3S1 = 


-Sl S3 


= s 2 


S 2 S 3 = 


-s 3 s 2 


= Si 


S 4 S k = 




= s k 



These rules are satisfied, in particular, by 
the Pauli spin matrices (only the first three 
bear this name, because 04 serves to form the 
identity matrix) [18], [20], [3] 





"0 


l" 




"0 


— i 


a 1 = 


1 





, 0"2 = 


i 







"1 


" 




i 


0" 


03 = 





-1 


, CT4 = 





i 



(11) 



where "1" in fflQj) represents the identity 
matrix, i the imaginary unit and Sj = — i o"j, 



for 16 {1,2,3,4}. 

The Pauli quaternions evidence one differ- 
ence with respect to the classical Hamilton 
quaternions, being the need of matrices, 
which in some cases have imaginary units i. 

The sum of two Pauli quaternions is of 
the same form as the given for the Hamil- 
ton quaternions and its product, using (fTOj) . 
becomes: 

A * * B =s^{ab — a x b x — a y b y — a z b z ) + 
Si(a6 x + a x b + a y b z - a z b y ) + 
s 2 {ab y - a x b z + a y b + a z b x ) + 
ss{ab z + a x b y - a y b x + a z b). 

(12) 

In a compact form, the product for the 
Pauli quaternions has exactly the same form 
as the Hamilton quaternions and, therefore, 
have the same problems as these: 

A**B = (ab— a-b, ab + b&+ axb) (13) 

1.4 Dirac matrices 

The Dirac matrices must satisfy the Klein- 
Gordon equation, the following relations 
should be satisfied by the Dirac matrices [55] : 

otiOij + atjOti = 25 ij: 

ai(3 + (3ai = 0, (*= 1,2,3) (14) 
a\ = ^ = I 

where / represents a N x N unit matrix. 
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For 2x2 matrices, only three anti- 
commuting matrices exist (the Pauli matri- 
ces) . Thus the smallest dimension allowed for 
the Dirac matrices is N = 4. If one matrix 
is diagonal, the others can not be diagonal or 
they would commute with the diagonal ma- 
trix. We can write a representation that is 
hermitian (a matrix is hermitian if it is equal 
to the conjugate of its transpose), traceless 
(trace equal zero), and has eigenvalues of ±1: 

«*=(° o)' (* =1 ' 2 ' 3 ) ( 15 ) 

and 

(16) 

where oi are the 2x2 Pauli matrices and 
/ is the 2x2 unit matrix. 

Finally, we are ready to define the Dirac 's 
gamma matrices out of «j and /3: 

7° = P, 7* = Mi = 1,2, 3) (17) 

These matrices satisfy the relations: 

( 7 °) 2 = 1, ( 7 T = -1, (18) 

and all four matrices anticommute among 
themselves. These relations are compara- 
ble to the Hamilton basis fl2D-flSD an< ^ P am i 
basis (flOl) . except for the exchange of sec- 
ondary importance in the signature of the 
gamma matrices, (+,—,—,—), to the sig- 
nature satisfied by the Hamilton and Pauli 
bases: (-,+,+,+)• 



2 The fourvectors 

The fourvectors are four- dimensional num- 
bers of the form 

A = e a t + i a x + j a y + k a z (19) 

or, assuming that the order of the basis el- 
ements is the indicated, those basis elements 
can be suppressed and included implicitly in 
a notation similar to a vector or 4D point: 

A = (a t ,a x ,a y ,a z ) (20) 

Where the elements of the fourvector can 
be any integer, real, imaginary or complex 
numbers. 

The four basis elements e, i, j, k satisfy the 
relations: 

e 2 = i 2 =j 2 = k 2 = e = eijk (21) 

The following rules are satisfied by the basis 
elements: 

e i = — i e = i, 

e j = -j e = j, 

e k = — k e = k, 

. . ' (22) 
1 J = -J i = k, 

j k = -k j = i, 
k i = -i k = j. 

The group of relations (1221) gives an im- 
portant operational mechanism to reduce any 
combination of two or more indices to at most 
one. 

The "e i j k" bases characterize the 
fourvector product as not commutative but, 
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what is more important and different with 
respect to the previous Hamilton and Pauli 
quaternions as well as to the Clifford Algebra 
(see [1], p. 5 axiom 3), the product is not 
always associative. For example consider 
the product of the following four symbols 
"iejk", in the order given; if, with the use 
of (I22l) . we first reduce "i e" to — i then 
"-i j" to -k and finally "-k k" to -e we 
obtain one result; but, if we first reduce 
the two middle basis elements "e j" to j, 
then "i j" to k and then "k k" to e, we 
get the same result but with the sign changed. 



If we put these rules into a multiplication 
table they look in the following way: 



** 


e 


i 


j 


k 


e 


e 


i 


j 


k 


i 


i 


e 


k 


-j 


j 


-j 


-k 


e 


i 


k 


-k 


j 


i 


e 



Let us assume that each basis unit is 
affixed its proper number (as components of 
a tensor) : 

(e,i,j,k) (90,91,92,93) = 9, 
then the fourvector multiplication satisfies 
the following compact relation, where the 
symbol \ has some similitude with the 
Levi-Civita symbol. Refer to [27] for the 
meaning of such symbol. 

QiQk = S ik q + xlklj h3> fc = 0, 1, 2, 3 



2.1 Discussion 

Tait and Gibbs discussed about the rela- 
tive merits of the "vector product" against 
the "quaternion product". Tait appeared to 
gain a slight advantage by pointing out that 
quaternion products are associative, whereas 
the cross product is not. Tait used the follow- 
ing example to reveal the supposed deficiency 
of the vectors [32] : 

i x (j x j) = ^ (i x j) x j = -i (23) 

Note that the four- vector product proposed 
in the previous section, although not associa- 
tive, or precisely because of that, gives the 
correct result for the problem at hand: 

i x (j x j) — > i x e —> -i / s 

(24) 

(i x j) x j -> k x j -> -i 

The fourvectors have extensive applica- 
tions in electrodynamics and relativity. The 
present author believes that the use of the 
fourvectors, with the proposed algebra, can 
replace advantageously the matrices, vec- 
tors and tensors in representation. Some of 
the advantages proposed for the Hamilton 
quaternions, Geometric Algebra and Space- 
Time Algebra, which are also extended to the 
fourvectors, are: 

1. Fourvectors can express rotation as a ro- 
tation angle about a rotation axis. This 
is a more natural way to perceive rota- 
tion than Euler angles [5]. 

2. Non singular representation (compared 
with Euler angles, for example) 
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3. More compact (and faster) than matri- 
ces. For computation with rotations, 
fourvectors offer the advantage of requir- 
ing only 4 numbers of storage, compared 
with 9 numbers for orthogonal matri- 
ces [24] . Composition of rotations re- 
quires 16 multiplications and 12 addi- 
tions in fourvector representation, but 27 
multiplications and 18 additions in ma- 
trix representation... The fourvector rep- 
resentation is more immune to accumu- 
lated computational error. [24]. 

4. Every fourvector formula is a propo- 
sition in spherical (sometimes degrad- 
ing to plane) trigonometry, and has the 
full advantage of the symmetry of the 
method [30] . 

5. Unit fourvectors can represent a rotation 
in 4D space. 

6. Fourvectors are important because of 
their "all-attitude" capability and nu- 
merical advantages in simulation and 
control [28] . 

Quaternions have been often used in com- 
puter graphics (and associated geometric 
analysis) to represent rotations and orienta- 
tions of objects in 3D space. This chores 
should be now undertaken by the fourvec- 
tors, which are more natural, and more com- 
pact than other representations such as ma- 
trices, and operations on them such as com- 
position can be computed more efficiently 
Fourvectors, as the previous quaternions, 
will see uses in control theory, signal pro- 
cessing, attitude control, physics, and or- 



bital mechanics, mainly for representing rota- 
tions/orientations in three dimensions. The 
spacecraft attitude-control systems should be 
commanded in terms of fourvectors, which 
should also used to telemeter their current 
attitude. The rationale is that combining 
many fourvectors transformations is more nu- 
merically stable than combining many matrix 
transformations. 

2.2 Complex fourvectors 

The only difference with respect to the 
ordinary fourvectors is that the elements are 
not purely real but complex numbers. 
The collection of all complex fourvectors 
forms a vector space of four complex 
dimensions or eight real dimensions. Com- 
bined with the operations of addition and 
multiplication, this collection forms a non- 
commutative and non-associative algebra. 
There is no difficulty in obtaining the mul- 
tiplicative inverse of a complex fourvector, 
when it exists, within the fourvector algebra 
suggested below. However, there are complex 
fourvectors different from zero whose norm 
is zero. Therefore the complex fourvectors 
do not constitute a division algebra. 

It is important to realize that the relations 
needed by the Klein-Gordon equation (fUj) . 
are directly satisfied by the purely real 
fourvectors, whereas the relations needed 
by the Dirac equation (TT8"j) , are satisfied 
by the fourvectors constituted of imaginary 
components in the vector part. 

This seems to suggest that there are two 
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different kinds of physical entities, although whose scalar component is equal to zero), 
closely related, which need respectively the Divided by two serves to define the corn- 
real and the imaginary representations. This mutator (similar to the vector Hamilton's 
insight appears potentially useful for Physics, operator V): (A — A)/2 = V A 



3 Fourvector algebra 

The sum of two fourvectors is another 
fourvector where each component has the 
sum of the corresponding argument compo- 
nents: 



A + B =e(a t + b t ) + i(a x + b x ) + 
j(oy + by) + k(a z + b z ) 



(25) 



The difference of two fourvectors is defined 
similarly: 



A - B =e(a t - b t ) + i(a x - b x ) + 
j(a y - by) + k(a z - b z ). 



(26) 



The conjugate of a fourvector changes the 
signs of the vector part: 



A = ea t 



}a y - ka z 



(27) 



From this definition it is obvious that 
the result of summing a fourvector with 
its conjugate is a fourvector with only 
the scalar component different from zero. 
Divided by two, isolates the scalar compo- 
nent of a fourvector and serves to define 
the operator named the anticommutator 
(similar to the scalar Hamilton's operator 
S): (A + A)/2 = SA. Similarly, the result 
of subtracting the conjugate of a fourvector 
from itself is a pure fourvector (that is, one 



The complex conjugate or hermitian conju- 
gate of a fourvector changes the signs of the 
imaginary parts. Given the complex fourvec- 
tor: 



A =e(a t + ib t ) + i(a x + ib x ) + 
j(a y + ib y ) + k(a z + ib y ) 
(10.15) 

Then its complex conjugate is: 
A 



f28) 



=e(a t - ib t ) 
j(a y -ib y ] 



\{a x - ib x ) + 
k(a z - iby) 



(29) 



Using relations (1211) and fl22|) . the fourvec- 
tor product is given by: 

A * * B =e(a t b t + a x b x + a y b y + a z b z ) + 
i (a t b x - a x b t + a y b z - a z b y ) + 
j (a t b y - a x b z - a y b t + a z b x ) + 
k(a t b z + a x b y - a y b x - a z b t ). 

(30) 

Using the notation of three-dimensional 
vector analysis we obtain a shorthand for the 
product. Regarding i, j, k as unit vectors in 
a Cartesian coordinate system, we interpret 
the fourvector A as comprising the scalar a 
and the vector part a = i & x + j & y + k a 2 . 
Then we write it in the simplified form A = 
(a, a). With this notation, the product (fBTii) 
is expressed in the compact form: 

A**B = (a 6 + a-b, ab-afr + axb) (31) 
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The following properties for the product 
are easily established: 

1. If the scalar terms of both argument 
fourvectors of the product are zero then 
the resulting fourvector contains the 
classical scalar and vector products in 
its respective components. 



2. The product is non- commutative. So, in 
general, there exist P and Q such that 
P**Q ^ Q i\ 

3. Fourvector multiplication is non- 
associative so, in general, P**(Q**R) 
^ (P**Q)**R 

Note that this is different from the 
Hamilton quaternions and the so-called 
Clifford Algebras, see for example pQ. 

4. The product of a fourvector by itself pro- 
duces a result different from zero only in 
the first or "scalar" component, which is 
identified as the norm of the fourvector. 
In this sense this constitutes the classical 
dot product of vector calculus: 

A**A = (a 2 + a 2 x + a 2 y + a 2 z , 0,0,0) (32) 

Note that this expression is substantially 
different with respect to the Hamilton 
quaternions, in which the square of a 
quaternion is given by 



(a 2 -v- v, 2a t v), 



(33) 



scalar component has terms with the 
sign changed, but appears a non-zero 
term in the vector part of the quaternion. 
This has been a source of difficulty to 
apply Hamilton quaternions in Physics, 
which is overcome by the fourvectors. 

5. The multiplicative inverse of a fourvector 
is simply the same fourvector divided by 
its norm. 

6. Properties of the product and conju- 
gates: 

P**Q = Q * * P (34) 



P**(Q**R)=R**(Q**P) (35) 

p * *(P * *Q) = P * *(P * *Q) 
= |P| *Q 

(P * *Q) * *P = (P * *Q) * *P 

= |P| *Q 



(36) 
(37) 



where v represents the three vector 
terms of the quaternion. Not only the 



With an operator notation: The prod- 
uct of two fourvectors is equal to the 
conjugate of the same product in reverse 
order: A**B = Conjugate[B**A] 



7. For the case that "r" is a rotor (a 
fourvector with |r| = 1) then: 

r * *(r * *Q) = Q = r**(r* *Q) 

(r * *Q) * *r = Q = (r * *Q) * *r 

((((Q * *r) * *r) * *r) * *r) = Q 

(otherwise, if |r| is not equal to 1, the 
products of this numeral are equal to Q 
or Q multiplied by |r|.) 



11 



8. The fourvectors do not contain the 
complex numbers, as is usually demon- 
strated for the Hamilton quaternions. 
The product of the fourvectors: (a, b, 0, 
0) and (c, d, 0, 0) is 

(ac+bd, ad-bc, 0, 0); also, the product 
of the fourvectors (a, i b, 0, 0) and (c, i 
d, 0, 0) is (ac-bd, % ad-bc, 0, 0), whereas 
the product as complex numbers should 
be: (ac-bd, i ad+bc, 0, 0) 

9. Given the fourvectors A and B, the 
commutator: 

[A, B] = i(A * *B — B * *A) (38) 
= (0, ah - ab + a x b) (39) 

gives a fourvector with zero scalar and 
with the vector part equal to the vec- 
tor part of the fourvector product A**B. 

For the curious, this commutator satis- 
fies the properties of antisymmetry and 
linearity. The Jacobi identity is satisfied 
only for pure fourvectors. 

10. Given two fourvectors, A and B, the 
anticommutator: 

<A,B> = i(A**B + B**A) (40) 
2 

= (a 6 + a-b,0,0,0) (41) 



A**B and with the vector part equal to 
zero. 

11. Product is left distributive over sum: 
a**(b + c) = a**b + a**c 

12. Product is right distributive over sum: 
(a + b)**c = a**c + b**c 

13. The product of three "pure" fourvectors 
(defined as those whose scalar compo- 
nent is zero) can be expressed with the 
following vector products: 

a * *(b * * c) = 

(a • (b x c), a x (b x c) - a * (b • c) ) 

Where "•" and "x" are the standard 
vector dot and cross products and "*" 
represents the product of the scalar 
(b • c) by the vector a. The scalar com- 
ponent of the result can be recognized 
as the volume of the parallelepiped 
having edges a, b and c. Consequently, 
if the three vectors a, b and c are in 
the same plane (or parallel to the same 
plane) then the scalar component of the 
resulting fourvector product is zero. 

14. The following identity is also satisfied: 
(a**b) ** (a**b) = (a**a) ** (b**b) 

3.1 Product with matrices 



gives a fourvector with the scalar equal Given two fourvectors, p and q: 
to the scalar of the fourvector product p=(p0,pl,p2,p3), 
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q=(q0,ql,q2,q3), 
their product can be obtained as 
R=p**q 

The same product can be obtained multiply- 
ing the following matrix P by the (four)vector 
q: 





Po 


Pi 


P2 


P3 
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R = 


-Pi 


Po 


~P3 


Pi 




Qi 




-P2 


P3 


Po 


~Pl 








-Pa 


-P2 


Pi 


Po 
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The P matrix has the property 
P • P T = P T • P = I, where T represents 
the transpose and I is the identity matrix. 
(More precisely the diagonal elements have 
the form: Po 2 +Pi 2 +P2 2 +P3 2 , which are equal 
to 1 only if p is a unit fourvector; else, in 
the diagonal is obtained the norm of the p 
fourvector) . 

3.2 The norm 

The norm of a fourvector is defined by 

\(a t ,a x ,a y ,a z )\ = a 2 + a 2 + a 2 + a 2 (42) 

It can be computed as the scalar component 
of the product of the fourvector by itself. 
The norm satisfies the properties 

|A| = |A| (43) 



|P * *Q| = |Q * *p| = |P| * |Q| (44) 

The last property allows to conclude the 
following form of the four-squares theorem: 



(a 2 + a 2 + a 2 , + a 2 3 )(b 2 + b\ + b\ + b 2 3 ) = 
(a b + ai bi + a 2 b 2 + a 3 b 3 ) 2 + 
(a h -a 1 b + a 2 b 3 - a 3 b 2 ) 2 + 
(o b 2 - ai 6 3 - a 2 b + a 3 b{) 2 + 
(a 63 + «i h - a 2 bi - a 3 b ) 2 

(45) 

3.3 Identity fourvector 

The identity fourvector, let us denote with 1, 
has the scalar part equal to 1 and the vector 
part equal to zero: 1 = (1, 0, 0, 0). 
It has the following properties, where "p" is 
any fourvector: 
l**p = p 

p**l = p 

As you can see, this is the left identity. It is 
possible to find the right identity of a fourvec- 
tor but it is a little more complex. See the 
section "Right factor of a fourvector" . 

3.4 Multiplicative inverse 

The multiplicative inverse or simply inverse 
of a fourvector P is denoted by P^ 1 . 

The inverse of a fourvector P is the same 
fourvector divided by its norm: 

P-^P/IPI (46) 

The inverse operation satisfies the proper- 
ties: 

p**p-i_-^ 
p— i**p_-^ 
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(p-l)-l = p 

(P ** Q)- 1 = P 1 ** Q" 1 
P _T = (P)- 1 
Commutativity of products including 
inverses: 

p ** (p-1 ** = p-1 ** (p ** Q) 

^p-1 ** ** p _ ** ** p-1 

Q = p**(p-l ** Q) = p~ T**(p ** Q) 

Q = (P- 1 ** Q)**P = (P ** Q)**^ 



3.5 Scalar multiplication 

Scalar multiplication If c is a scalar, or a 

scalar fourvector, and q=(a, v) a fourvector, 

then cq = (c, 0) * *q = (c, 0) * *(a, v) = 

(ca + • v, cv — aO + x v) 

Simplifying: 

c (a, v) = (ca, cv) 

3.6 Unit fourvector 

A unit fourvector has the norm equal to 1. 
It is obtained dividing the original fourvector 
by its magnitude or absolute value, that is 
the square root of the norm. The product of 
two unit fourvectors is a unit fourvector. A 
unit fourvector can be represented with the 
use of trigonometric functions 

w = (± cos(a) ± u sin(a)) 

where u is in general a 3D vector of unit 
length. 

The product of two unit fourvectors: 
Assume the unit fourvectors a and b: 



a = (cos(a), sin(o;), 0, 0) 
b = (cos(/3), sin(/3), 0, 0) 

Its product is 

a**b = (cos(/3 - a), sin(/3 - a), 0, 0) 
so, if a = b, then the resulting fourvector 
is the identity fourvector. 

The inverse of a unit fourvector is the same 
unit fourvector. This is because the product 
of the fourvector by itself produces the iden- 
tity fourvector, or the norm, equal to 1, in 
the scalar component. 

3.7 Fourvector division 

The fourvector division is performed by 
multiplying the "numerator fourvector", P, 
by the "denominator fourvector", Q, divided 
by its norm (or rather multiplying P by the 
inverse of Q): 

P * *Q/|Q| = P * *Q 1 

If P and Q are parallel "vectors" (pure 
fourvectors or with scalar part equal to zero), 
then the division produces, in the scalar part, 
the proportion between both vectors. For ex- 
ample: 

P=(l, 2, 3, 4) 
Q=(3, 6, 9, 12) 

P**Q/|Q| = (1/3,0,0,0) 

3.8 Right factor of a fourvector 

Let us try to solve the following equation for 
"X" 

A == B ** X 
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Let us assume that A and B are constant 
fourvectors and X is an unknown fourvector 

A=(a0,al,a2,a3) 
B=(b0,bl,b2,b3) 
X=(x0,xl,x2,x3) 
Where the components can be integer, ra- 
tional, real or complex. 

The solution for X can be obtained with 
the expression 

X = B~ 1 **A (47) 

For example, let us try to determine what 
the value of X should be in order to satisfy 
the equation A == B ** X, when 

A=(7,l,-3,5) 

B=(l,3+t 5,2,-1) 

According to equation (T4T1) the solution is 

X = ( =2 _ il _9 -4 12 _ -53 _9_ _ -21 \ 

V 10 5' 10 5' 25 50' 25 50^ 

Replacing this solution into the product 
B ** X it can be verified that it reproduces 
the A fourvector. 



3.9 Left factor of a fourvector 

In a similar form, let us try to solve the fol- 
lowing equation, where the unknown "Y" is 
now a left factor of the constant fourvector B: 

A == Y ** B 



Y is obtained with the expression 

Y = Hprod[A, B" 1 ] (48) 

Where "Hprod" represents the product for 
the classical Hamilton quaternions or Grass- 
man product. 

For example, for the same fourvectors A 
and B from the previous example, let us 
apply the equation (j4"8l and find 

y — (± — j± 39 _ • 29 _ 19 , 11 11 ■ ,-21\ 

^25 50' 50 25' 25 ~ 50 ' 50 ^ 25' 

Replacing this solution into the product 
Y ** B it can be verified that it reproduces 
the A fourvector. 

Note that the left and right factors of some 
fourvector such as B are different, although 
both factors have the same norm and satisfy 
the equality: 

X* * Y" 1 == X" 1 * * Y 

Both products, B ** X and Y ** B are 
equal to A. 

Formula (jlSjl can be used to determine the 
single rotation fourvector that produces the 
same effect as a rotor. However, the results 
obtained are, in general, more complex than 
a classical rotor. Nevertheless, if we need 
the fourvector L, which produces the same 
rotation of the fourvector p as the rotor 
fourvector r, that is: 

L * * p = Rotate[p, r] = r * * (p * * r) 

applying equation ( 1481) we solve for L: 
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L = Hprod[(p * * r) * * r], p x ] 

Where "Hprod" is the Grassman product. 



3.10 Solution of quadratic 
fourvector polynomials 

There can be an infinite number of solutions 
for a quadratic fourvector polynomial. Con- 
sider the quadratic equation 

q 2 == 1. 

Then, the fourvector q = (cos(x), 
sin(x),0,0), where x is any real number is 
a collection of solutions for this equation 
because the norm of q is 1: 

q ** q == 1 

So the above choice for q satisfies the 
equation q 2 == 1 for all real values of x. 

When there is a solution of a quadratic 
equation, it can be computed as in the 
following. 

Assume a quadratic polynomial of the form: 
q**q+q**j==k 
where: 

q=(q0,ql,q2,q3) is an unknown fourvector 
and 

j=(0,-l,l,0) and k=(-l, 0,0,1) are constant 



fourvectors 

From here, the four equations, obtained 
equating the four components, are: 

1+qO 2 - ql+ ql 2 + q2+ q2 2 + q3 2 ==0, 

- q 0+q3==0, 

q 0+q3==0, 

-l-ql-q2==0. 

This system of equations has two solutions 
for the four components of the fourvector: 



q 1 = (0 ) -l-i/^ ) i/v^ I 0) (49) 
q 2 = (0, 1 + ^/v / 2, -i/V2,0) (50) 



Replacing q x (or q 2 ) by its value, in the 
following expression, which is the left hand 
side of the given equation, q**q+q**j, the 
value returned is: (-1,0,0,1), 

Which is the value of the right hand side. 

For a comparable quadratic equation, but 
now affecting the "j" fourvector to the left 
of "q" instead of to the right, 

q**q + j**q==k 

The two solutions are: 

qi = (0,-Vv / 2,l/2(-2 + zv / 2),0) (51) 
q 2 = (0,Vv / 2,l/2(-2-. v / 2),0) (52) 

Which, replacing in q**q+j**q 
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produce the value: (-1,0,0,1) 

Note that the fourvectors form a division 
algebra since they have a (left) multiplicative 
identity element 1^0 and every non-zero 
element a has a multiplicative inverse (i.e. 
an element x with ax = xa = 1). 



4 Fourvector rotation 

Mathematically, a rotation is a linear trans- 
formation that leaves the norm invariant. 
These are called orthogonal transformations. 

There are several methods to represent 
rotation, besides quaternions, including 
Euler angles, orthonormal matrices, Pauli 
spin matrices, Cayley-Klein parameters, and 
extended Rodrigues parameters. 

From the several ways to represent the 
attitude of a rigid body one of the most 
popular is a set of three Euler angles. 
Some sets of Euler angles are so widely 
used that they have names, such as the 
roll, pitch, and yaw of an airplane. The 
main disadvantages of Euler angles are 
that certain functions of Euler angles have 
ambiguity or singularity for certain angles. 
This produces, for example, the so-called 
"gimbal lock". Also, they are less accurate 
than unit quaternions when used to integrate 
incremental changes in attitude over time [7]. 

The handling of rotations by means of 
quaternions has constituted the technical 



foundation of modern inertial guidance sys- 
tems in the aerospace industry for the orien- 
tation or "attitude" of satellites and aircrafts. 
This task is to be shown in the following that 
should be carried out by fourvectors. 

Many graphics applications that need to 
carry out or interpolate the rotations of 
objects in computer animations have also 
used quaternions because they avoid the 
difficulties incurred when Euler angles are 
used. The form to replace by fourvectors is 
not performed in the present paper, although 
it should be perfectly possible. 

The use of matrices is neither intuitive for 
the localization of the axis of rotation nor 
efficient for computation. But one of the 
most important disadvantages is the asso- 
ciativity of both the matrix and Hamilton's 
quaternion products where, for example, 
A • (B ■ C) is equal to (A • B) • C. In fact it 
is rather well known that the composition of 
rotations, when either matrices or Hamilton 
quaternions are used, is associative. This 
means that these mathematical tools produce 
the same result no matter what the grouping 
of a sequence of two or more rotations. This 
poses a serious technical problem to the 
engineers who need to distinguish between 
two sequences of the same rotations. To 
illustrate with an example, let us assume 
that you are piloting an airplane with a local 
frame of reference whose origin is attached 
at the center of the airplane. Assume that, 
at a certain instant, the "x" axis, which is 
pointing toward the front of the plane, is 
horizontal according to an observer in the 
earth, let us say directed toward the North 
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pole, the "y" axis points to the right wing, 
that is pointing to the East, and the "z" axis 
points toward the earth's center. If under 
these conditions you maneuver to produce a 
90° roll (a quarter circle clockwise rotation 
about the local x axis) followed by a 90° 
yaw (rotation "to the right" around the 
local z axis) your plane should be falling 
perpendicularly to the earth; but if you 
interchange the order of these rotations, 
that is first the yaw and then the roll, 
this should put your plane in a horizontal 
heading toward the East, although with the 
right wing pointing toward the earth's center. 

To rotate a vector u around an Euler vector 
n through an angle 9, one has to apply the 
following equation (see Fig. 1 and references 

[si, my 

v = ucos9 + n(n • u)(l — cos9) + (n x u)sin9 

(53) 



o 




Figure 1: Graphic of a rotation 

According to Silberstein, [26]: "It has been 
remarked by Cayley, as early as in 1854, that 



the rotations in a four-dimensional space may 
be effected by means of a pair of quaternions 
applied, one as a prefactor and the other as 
a postfactor, to the quaternion u whose com- 
ponents are the four coordinates of a space- 
point, say v = aub". This phrase applies 
directly to fourvectors if "quaternion(s)" is 
replaced by "fourvector(s)". 
In the case of pure rotation, a and b must 
be either unit-fourvectors or the norm of their 
product must be 1: |a| * |b| = 1. 
It follows, from the rule: 

|a * * b| = |a| * |b|, (54) 

the multiplication of the fourvector being ro- 
tated by unitary fourvectors a and b, effects 
an orthogonal transformation. 
This form can be simplified so instead of two 
different unitary fourvectors is selected only 
one, let us name it r. 

A possible fourvector r that produces the 
rotation of any fourvector V about a certain 
axis "n", through an angle 9, has the form 

ra, 0, m na, usi, my 

r = (cos(^),n x sin(-),n y sin( 6 -),n z sin(-)) 

(55) 

or simply 

9 9 

r = (cos(-),n sin(~)) (56) 

The rotation is carried out with the following 
product: 

V' = r**(V**r) (57) 

The rotation operand needs to be conjugated. 
This is different with respect to the rotation 
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using Hamilton quaternions [2B], where the 
inverse or the conjugate of the second rotor 
is needed. 

Fig. 2 shows the vector V rotated around 
the unit vector k, through an angle 8, with 
rotor r. The rotation operator is linear. It 



q . * * q t 

x axis: 



1, and cause rotations about the 




Figure 2: Example rotation 



can be proved that the rotation of either the 
product or the sum of two fourvectors s and t 
with the rotor r, is the same as, respectively, 
the product or the sum of the rotations of 
each fourvector: 

Rotate [s * * t, r] = Rotate[s, r] * * Rotate[t, r] 

(58) 

Rotate[s + t,r] = Rotate [s,r] + Rotate[t,r] 

(59) 

Let us define the following example rotation 
fourvectors, or rotors, whose norm is the unit: 



qi = (±cos(a/2), ±sin{a/2), 0, 0) (60) 
q 2 = {±cosh(a/2), ±i sinh{a/2), 0, 0) (61) 



03 



7 + 1 . / 



,0,0) 



(62) 



2 ' V 2 
where 7 is the Lorentz contraction factor. 

The rotation fourvector q 3 was obtained by 
transforming the following fourvector 
(±7, ±i/3 7 , 0, 0) 

in such a way that it produces a rotation 
of half the hyperbolic angle, as with q 2 . 

Let us multiply any one of the previous 
fourvectors by the following one that repre- 
sents a differential of interval: 

ds = (c dt, i dx, i dy, i dz) 

The products are of the form: 
Rotate[ds, qj = q^ * * (ds * * q i ) (63) 

where q« is anyone of the above list of unit 
fourvectors. 

Then, the norm of the result is the square of 
the differential of interval: 

ds 2 = c 2 dt 2 - dx 2 - dy 2 - dz 2 

This means that any one of these trans- 
formations (rotations) preserves the interval 
invariant. 

But first the rotation of a real fourvector 
a = (a t ,a x ,a y ,a z ) with the rotor q 1 produces 
the classical formulas for rotation of a vector 
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about the x axis: 

(a t ,a x cos(2a) - a y sin(2a), 
Rotate^, Ql ] = (a t , a x , ^ cog(2a) + ^ ^ (67) 



a y cos(a) — a z sin(a), (64) 
a z cos(a) + a y sin(o;)) 



4.1 Composition of rotations 



The opposite rotation can be done using The rotation through an angle a followed 

as rotor the inverse of q 1? which changes by another rota ti n through an angle (3 is 

the sign of the angle a. For example, if we equiva i ent to a single rotation through an 

rotate the last result with the conjugate or ari gi e a + ft- 
the inverse of q 1 then the original fourvector 
a is recovered. 



The rotation of the differential of interval 
with q 2 gives the complex fourvector: 



Assume, for example that the rotations are 
produced by application of the following ro- 
tors: 



rothl = (cosh(a/2),i sinh(a/2), 0, 0) (68) 
roth2 = (cosh(P/2),i sinh(3/2), 0, 0) (69) 



Rotate [ds, q 2 ] = (cdt, idx, 

dz sinh(a) + idy cosh(o;), (65) 

— dy smh(a) + idz cosh(a)) L e t, us apply these rotors over the fourvec- 

tor M=(a,b,c,d) with the operation: 

The rotations with q 3 produce Lorentz 
boosts. Let us apply to ds: M1 = Rotate ^ roth ^ ( 70 ) 

Rotate [ds, qJ = (cdt, idx, c ,, , , f u ■ , ,. 

L ' ^ J v ' ' followed by the following rotation: 

i 7 (dy — i/3 dz), (66) 

17 (dz + ifidy)) M2 = Rotate[Ml, roth2] (71) 



If the fourvector to be rotated is the 
previous a and the rotor fourvector is 



We obtain the following result: 



(a, i b,i ccoshict + 3) + d sinh(a + 3), , 
r = (cos(a/2),0,0,sin(a/2)) - A u) Js -u) Js ( 72 ) 



idcosh(a + 3) — csinh(a + 3)) 

Which is identical to the rotation produced 
by directly applying over M the rotor: 



Then the double rotation: 
r * *(r * *(a * *r) * *r) 

Is equal to a single rotation through the roth = [cosn( - ),tstnh{ - ),U,UJ 
double angle 2a. The result is: (73) 
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4.2 Reflections 

Let us show how the reflection in a plane 
with unit normal a can be done (see next 
figure) 




Figure 3: Reflection 



the normal vector a, which defines the 
direction of the plane, generates a reflection 
of the vector x if we rotate x around the 
vector a through an angle of 7r radians and 
then change every sign of x (otherwise we 
end up with an arrow pointing at the same 
point as x). This rotation is done by fixing 
at 7r radians the angle a of rotation of a 
rotor of the form ql of section HJ i.e. 

q[Cos[a/2],a Sin[a/2]]. 

Consequently, the cosine term disappears 
and the sine term becomes equal to 1, with 
which we are left with the vector a alone, as 
rotor. 

The vector x, after reflection, is: 
x' = —a * *(x * *a) 

To simplify this expression, let us note that 
the rotation through it radians clockwise is 
the same as the rotation through tt radians 



counterclockwise, so we could include conju- 
gations of the vector a. Also, the conjugation 
of the vector x can be canceled with the neg- 
ative sign, so the final equation is: 
x' = a**(x**a). 

As was suggested at the end of the previ- 
ous section, the rotations can be composed. 
So that, by multiple application of this re- 
flection formula it is possible to compute the 
vector x' that describes the direction of emer- 
gence of a ray of light initially propagating 
with direction x and reflecting off a sequence 
of plane surfaces with unit normals ai; 3l2] 




Figure 4: Reflections 



Hestenes [15] shows other applications 
for rotations and reflections, for example 
in crystals. Note that in his "geometric 
algebra" Hestenes needs the negative sign 
that we don't need. 
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5 Conclusions 



The mathematical structure of quaternions 
has always been considered as more appro- 
priate than the simple vectors to represent 
the real physical variables. Nevertheless, 
the quaternions were dismissed for the 
difficulties and complications produced by 
their quaternion product. With the new 
product, suggested in the present paper for 
fourvectors, all those difficulties disappear. 
Of course there must be a delicate balance 
between the correct mathematical tools and 
the real physical objects being studied and 
handled. One has to also be aware that 
mathematics clearly affects the ontology of 
physics [TT] . 



The fourvector algebra proposed in the 
present paper seems to be the correct 
mathematical tool to study the fundamental 
physical variables and their describing equa- 
tions. 

This new mathematical structure is an 
extension of the classical vectors. Its sim- 
plicity contributes to the possibility of more 
extended and fruitful use in all branches of 
science. 

The applications of the Hamilton quater- 
nions for rotations in three dimensions 
have been the more extended in current 
Physics. Such uses, as well as reflections, 
are still permitted by the fourvectors. The 
applications to Lorentz boosts had problems 
with the old quaternions, so this area opens 
up for the scientists. 
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